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Abstract 

New data on polarized ^i—p and e—p scattering permit a first determination of a^ 

using the Bjorken sum rule, as well as liigher precision in determining the nucleon spin 

, decomposition. Using perturbative QCD calculations to ©((as/vr)^) for the non-singlet 

Q \ combination of structure functions, we find as(2.5 GeV^) = 0.375lQQg]^ , corresponding to 

1-^ ' UsiM"^) = 0.122^0009 ' ^^^ using calculations to 0((as/7r)^) for the singlet combination 

we find Au = 0.83 ±0.03, Ad = -0.43 ±0.03, As = -0.10 ±0.03, AS = Au + A(i± As = 

^ . 0.31 ±0.07, at a renormalization scale Q^ = 10 GeV^. Perturbative QCD corrections play 

H \ an essential role in reconciling the interpretations of data taken using different targets. 

We discuss higher-twist uncertainties in these determinations. The Af; determinations are 

used to update predictions for the couplings of massive Cold Dark Matter particles and 

axions to nucleons. 
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When Bjorken first derived his sum rule for polarized deep-inelastic scattering |jl|, 
he doubted whether it could be tested experimentally. Several years later, however, 
he changed his mind |^ and an experimental programme on polarized electron- 
proton scattering was launched at SLAC 0,^,0 In view of this, a sum rule for the 
polarized proton and neutron structure function gi was proposed 0, based on the 
dynamical hypothesis that strange quarks and antiquarks in a polarized nucleon 
would have no net polarization. However, this sum rule was not at all rigorous, 
in contrast to the Bjorken sum rule which is an inescapable prediction of QCD 
[0,H. The first round of polarized e-p experiments that were consistent, within 
their stated errors, with the hypothesis of no net strange polarization. However, no 
measurements were made on polarized e-n scattering, so the crucial Bjorken sum 
rule went untested. 

Several years later, following the advent of naturally-polarized high-energy and 
-intensity muon beams at CERN, the EMC made a second measurement of Qi PJTO 



This was more precise than the earlier SLAC measurements, and extended to lower 
values of a; = \q^\/2mpV. It indicated a disagreement with the polarized-proton sum 
rule, corresponding to a non-zero and negative net contribution As of strange quarks 
and antiquarks to the spin of the proton, and a small total contribution of the light 



quarks: An + /S.d -\- As <^ 1 [in|'0- This result was surprising from the point of 
view of the most naive formulation of the constituent quark model, which would 
suggest Am + Arf ^ 1, and As = 0, and even compared with more sophisticated 
versions adapted to fit measurements oi qa = Am— Ac? and the hyperon axial-current 
matrix elements' F/D ratio, which suggested Am -|- Ad ~ 0.6 and As ~ 0. 

This surprise indicated that our theoretical understanding of non-perturbative 
QCD was incomplete, and stimulated attempts to remedy this defect. For example, 
it was pointed out |jl3| that the EMC result could be understood qualitatively within 



SU{3)f topological models of the nucleon, which describe light quarks via collective 
fields, incorporate the global symmetries of non-perturbative QCD, and predict Am+ 
Ad+As <^ 1. Alternatively, it was suggested [0,|l^,|l^ that the f/(l) axial anomaly 
might play a key role, which would modify the naive quark model predictions. It 
was even suggested that the "sacred" Bjorken sum rule might be violated P 



In view of the interest in checking the EMC polarized fi-p results and extending 
them to test the Bjorken sum rule, extensive experimental programmes are under 
way at CERN, SLAC and DESY. An important round of results on polarized ii- 
D scattering from CERN |jl8| and polarized ^He scattering from SLAC [|l^ were 



published in 1993, permitting for the first time an experimental test of the Bjorken 
sum rule. The CERN and SLAC data agreed within their errors on the extracted 
value of Qi and confirmed the validity of the Bjorken sum rule with a precision of 



about 12% ||2^, once the then-available higher-order perturbative QCD corrections 



np were taken into account and allowance made for higher- twist corrections |2^J2^ . 
These data also indicated As < and Am + Ad + As ^ 1. 



Recently, two new sets of data with polarized proton targets have been made 
available [0,^. In addition, order (a^/vr)^ corrections to the singlet part of struc- 
ture functions gf and Qi have recently been calculated |^ and estimates made of 
higher-order perturbative QCD corrections to the Bjorken and singlet sum rules 
2^,^. Therefore, now is an appropriate moment to reassess the precision with 



which the Bjorken sum rule has been tested, and to extract new estimates of Am, 

Ad and As. 

We find that the Bjorken sum rule is now verified with a precision of 10%, 
once all available perturbative corrections are taken into account. As an exercise, 
we determine a^ for the first time from the polarization sum rule data, finding 
a,(2.5 GeV^) = 0.375^^8? , corresponding to a^iMl;) = 0.122t[|;gE]^ . On the other 
hand, we find a consistent pattern of violation of the separate sum rules for proton 
and neutron targets 0, indicating that As < 0. Encouraged by the consistency, we 
extract values of Am + Ad + As from the different data sets. These are superficially 
different if the data are analyzed in the naive parton model, i.e. neglecting perturba- 
tive QCD corrections. However, the agreement between different data sets improves 
systematically as each higher order of perturbative QCD correction is included. 
Including the C((as/7r)^) calculation |21|, together with estimates of 0{{as/n)^) 
effects in the non-singlet channel ||2^ and the 0{{as/n)^) calculation |^, together 
with estimates of (!?((as/7r)^) effects in the singlet channel [^], we find 

An = 0.83 ±0.03, Arf = -0.43 ± 0.03, As = -0.10 ± 0.03, AS = 0.31 ± 0.07 

(1) 
at a renormalization scale Q^ = 10 GeV , with a global x^ = 1-3 ^oi 5 degrees 
of freedom. As an example of relevance of these results, we present at the end of 
this paper an analysis of the implications of these determination for dark matter 
couplings to matter, including the elastic scattering of supersymmetric relics and 
axion couplings to nucleons. 

Up to now, the prevailing attitude has been to use polarized structure function 
data to test the Bjorken sum rule, using the most complete theoretical calculations 



of corrections and a world-average value of a^. The conclusion last year was [^ 
that the EMC/SMC and SLAC E142 data together verified the Bjorken sum rule 
to within the available precision of 12%. Here we take a different attitude, more 
akin to that adopted with regard to the Gross-Llewellyn Smith sum rule |^^. There 
is no convincing theoretical evidence against the validity of the Bjorken sum rule, 
which is a solid prediction of QCD. Therefore, its validity can be assumed. Instead, 
one can use the new polarized structure function data to extract a value of ag^Q"^), 
whose consistency with other measurements is an a posteriori check on the validity 
of the Bjorken sum rule. 

As is well known, data at low Q^ are particularly sensitive to Ajg^ and hence 
at a premium in determining as(M|)| . As is illustrated by r decays, even a 

MS 

relatively imprecise determination of ^^(low Q'^)jjs (provided higher-order pertur- 



bative QCD and nonperturbative uncertainties are understood) extrapolates to a 
relatively precise determination of as{M'^)\ pl|. The SLAC E142 ^He (i.e. n) 

MS 

data have {Q"^) ~ 2 GeV^, and the E143 p data have (Q^) ~ 3 GeV^, and are hence 
particularly well placed to exploit this lever arm on Ajj^. Since they are at higher 
Q^, the new SMC p data are less useful in this respect, though they do provide 
important information about gi{x, Q"^) at small x, and enter into the determination 
of the Aq, as we discuss later. 



We use the value oiT'^{Q^=2 GeV^) = -0.028±0.006±0.009 given in ref. ^ and 
the value r?(Q2=3 GeV^) = 0.133 ± 0.004 ± 0.012 given in ref. ^. We emphasize 



that these estimates are based on polarization asymmetry Ai{x,Q ) data taken at 

average values of Q^ that depend on the bin in x, which have been converted into 

values of gi{x, (Q^)) by assuming that the Q^ dependence of Ai{x, Q^) is insignificant 

2^,0 and using standard parametrizations of F2{x,Q'^) |Q and Ri^x^Q"^) [^ to 



estimate 



gi{x,Q^ 



Ai{x)F2{x,Q^ 



(2) 



2x[l + R{x,Q^)] 

No data set indicates a significant Q^ dependence of y4i(x, Q^), and perturbative 
QCD calculations p5| , p6| , p7| , p8[ lead one to expect that the Q^ dependence can be 
neglected at the level of precision required here. The estimates of gi'"'{x, Q"^) also 
require assumptions on g2{x, Q"^) that are borne out by the latest experimental 

m 



bounds 



and consistent with the latest theoretical calculations 39 . Based on 



the above numbers, we use in our subsequent analysis 



pp_pn| 



=2.5 GeV^ 



0.161 ±0.007± 0.015 



(3) 



for the Bjorken integral. 

The small-x behaviours of gi'^ have recently been discussed ^^ in the context of 
a non-perturbative model of the Pomeron [^, which suggests that gf'"' ^ log(l/a;) 
at small x, rather than ~ x", —0.5 < a < |42,43|. Although such an effect 
would alter the estimates of F^''" that we use, it would cancel out in the difference 
that appears in the Bjorken sum rule, since this effect would be due to two-gluon 
exchange. Hence the estimate in equation (^ would be unaltered. 

Perturbative QCD corrections to the Bjorken sum rule have been calculated up 

to 0{{as/TT)^) |2T]] and an estimate was made of the C((as/7r)^) coefficient [ P^ : 



1 



T^^m-^iiQl = 7:\9A\ 
o 



l_r^^U 3.5833 f^ 



TT 



TT 



(4) 



_20.2153 I ^^^V - 0(130) f^^^ I + 



TT 



n 



in the MS prescription for Nf = 3 flavours, as appropriate to Q^ = 2.5 GeV^. Tlie 
value of as{2.5 GeV^) extracted by comparing (|^) and (^) depends on the order in 
QCD perturbation theory which is used. As we see in fig. 1, the extracted value 
decreases as one progresses from 0{as/7T) to higher orders, but shows good signs of 
stabilizing in (9((as/7r)^). We infer from this analysis a value 



a. 



[2.5 GeV')\jrs,N,=3 = 0-^75t'oZ (5) 



which corresponds to 



«.(M|)U,^^=5 = 0.122tno9 (6) 

Note that in evaluating (Mj we have used the 3-loop renormalization group equa- 



tions ^^, and matched values of a^ at the flavour thresholds Q = mc,mi„ as 



is appropriate in the MS scheme. The relevant prescription is given in ref. 
as-{mQ) = as^{mQ)[l + a{as+{mQ)/TT) + 6 (^^^(mQ)/??)^], where a = 0, b = 7/72. 
This means that the two-loop as-{mQ) = as_^_{mQ), but that the three-loop as-{mQ) 
and as+{mQ) are slightly different at mq. Figure 3 compares this determination of 
as with others, as compiled in p6|. 



We have not yet included higher-twist effects [^] in this analysis, for two main 
reasons. One is that their coefficients appear to be rather smaller than originally 
thought p3l, and also because the available estimates differ considerably from one 



another |47.48.4^,50,51,52,53|. For our purposes it is sufficient to take a rough 
estimate, 

Cht (-0.02 ± 0.01) GeV^ 

Q2 - g2 



S,, (F? - F^) = ^ ^ ^ w...^.^..;^.v ^^^ 



Another reason is that it may be double counting |Q to include this as well as the 



higher-order perturbative QCD contributions that have now been calculated and 
estimated. Including the estimate (0) in a fit to the data described above, we find 

«.(^l)lMS,iv,=5,H. = 0.118t°:°;i (8) 

The difference between this and the value in equation @ is indicative of the theo- 
retical systematic error. 

This Bjorken sum rule determination of as is encouragingly precise, and quite 
consistent with the other determinations shown in fig. 3, though not yet competitive 
with the market leaders. This consistency means that the Bjorken sum rule is verified 
within 10% precision by the SLAG E142 and E143 data alone, as well as being 



verified to within 12% precision by the SLAG E142 and EMG/SMG data 0. The 
higher precision expected from future SLAG data might enable this determination of 
as to become truly competitive. A good understanding of the radiative corrections 



will be needed at this level [^Sf. The optimal evaluation of the Bjorken integral will 



also require data from low x, where the SMG data are available and will continue 



to provide crucial input. Also, the very high precision data from HERMES |^ at 



moderate x will be very useful. The Bjorken sum rule is making the transition from 
a test of QCD to a tool for evaluating ag within QCD. 

Having verified that the Bjorken sum rule is well satisfied by the latest data, 
we consider the sum rules for proton and neutron targets separately. Their theo- 
retical prediction requires an extra assumption on the singlet axial current matrix 
elements AE((5^) or, equivalently, the nucleon matrix element As of the 37^755 
current. Alternatively, one can regard measurements of r^^"" = /q dxgi {x,Q^) 
as determinations of AZ1((5^), or equivalently As, if one assumes the validity of the 
Bjorken sum rule. This is the attitude taken here. The fact that Ff — V^ is in 
good agreement with the Bjorken sum rule means that the different targets yield 
consistent values of AS, as we now show. 

We use in our determinations of AE(Q^) and As the calculated 0{{as/TT)^) cor- 



rections to the Bjorken sum rule and the estimate [^ of the 0{{as/TT)'^) corrections 
used above. We also use the recent calculation of the 0((as/7r)^) correction [|^ 



to the singlet part of the sum rule and a recent estimate |^ of the 0{{as/7i)^) 
correction to it. If one chooses the renormalization scale /i = Q, the proton and 
neutron sum rules including these corrections can be written as 



X 




dxgy {x, Q 



3.5833 




1 , , 1 

± — \Qa\ H Os 

12'^ ' 36 



-20.2153 



'as{Q' 



■K 



-0(130) 



'a.(g' 



| + ... 



7r 



+ 



'asiQ' 



1.0959 



'«.(Q' 



TT 



TT 



0{6) 




-AS(Q^ 



(9) 



It is worth emphasizing that the Q^ dependence due to QCD in the singlet channel 
has two sources, namely the anomalous dimension of the singlet axial current, as 
well as the coefficient function. However, the Q^ dependence of AS((5^) is not large 
from the Q^ range of the present data upwards. Note that we have not included any 
higher-twist contributions, but will comment later on their possible effects. 

It is interesting to note that, because of their different Q^-dependences, it would 
be possible in principle to separate the SU{3) octet and singlet contributions (og, AS) 
to either Ti^Q"^) or r"((5^), and avoid in this way the use of SU{3) ||5^ relations 
to estimate as on the basis of hyperon /3-decay data. Writing the square-bracket 
perturbative QCD correction factors on the first and last lines of equation (^ as 
f{as) and h{as), respectively, for the octet and singlet contributions, one can write 



36 



+ 



h{as 



-AS(Q^ 



rW 



12 



-9A 



/M„, + M^AE(g=) = rjw^) + i^s, (10) 

l^a, + M^AE(e^) = rf(Q=) 



where the right-hand sides are combinations of measurable and theoretically-known 
quantities. (Strictly speaking, Ff in eq. (|TOD is (F^ + F")/2, i.e. it includes nuclear 
corrections.) In leading order, /(a^) = h{as) = 1 — (as/vr), and the different 
equations (JT^) are not independent. At NLO, however, /(a^) 7^ h{as), and then 
each of the equations (0) allows in principle an independent determination of AE 
and Og, provided data with sufficiently high precision are available at different values 
of Q^. One can also combine data from experiments with different targets. The 
current precision of the world data does not permit a meaningful separation of ag 
and AE using this technique. 

Using the known variation of asiQ"^) and assuming AE ~ 0.3 in eq. (0), one 
can estimate the experimental precision required to make this approach practical. 
For a range of Q^ between Q^=2 GeV^ and (5^=15 GeV^ we estimate that the right- 
hand side of eq. (|TUp would have to be known with a precision much better than 
±0.002, including both statistical and systematic errors. Currently E143 quotes [^ 
error values of ±0.004 (stat.) and ±0.012 (syst.), so the required improvement in 
precision would have to be very substantial, especially in the systematic error, unless 
for some reason the systematic error is approximately independent of Q^, in which 
case it would largely cancel out when comparing the r.h.s. of eq. (|10|) for the various 



values of Q . Pending future improvements in the experimental precision, for the 



rest of this paper we assume the usual SU{3) value ag = 0.601 ± 0.038 |2^, and we 
return now to the discussion of the existing data. 

In order to combine all the data available from different targets, we re-express 
the measurement on Deuteron and ^He targets in terms of the corresponding values 
of r^{Q^), assuming the Bjorken sum rule and incorporating all the non-singlet 
perturbative corrections /(a^) on the first line of equation (|^). The resulting values 
of TilQ"^) inferred from the E142 ^He and SMC D data are shown together with the 
EMC, SMC and preliminary E143 p data in fig. 4. We see no signs of convergence 
towards the naively-suggested value holding if As = 0. This situation contrasts 
with what has been found previously with the Gross - Llewelyn Smith and Bjorken 
sum rules, namely good agreement once the available perturbative QCD corrections 
are included. Also shown in fig. 4, to guide the eye, is the prediction for Ti{Q^) 
that would be obtained if As = —0.10 ± 0.04 which is highly consistent with the 
data. The inclusion of all available higher-order perturbative QCD corrections is 
important for this consistent picture to emerge. We plot in fig. 5 the values of AE 
that would be extracted from each experiment if one restricted one's analysis to 
include only low orders of perturbative QCD. The values of AE extracted using the 
naive parton model, i.e. (9((as/7r)°), are in poor agreement, particularly the E142 



neutron measurement done on ^He. However, the agreement improves significantly 
as one proceeds to Oi^asln), C((«s/7r)^)and (9((as/7r)^). (This last analysis includes 
the estimate of the 0{{as/'T^Y) corrections to the Bjorken sum rule [^ and of 



the 0{{as/Ti) ) correction [28| to the singlet sum rules.) The overall x of the 



global fit decreases systematically as each order of perturbative QCD is included, 
X^ = 12 (naiVe) ^ 4.2 {0{ajTi)) ^ 2.4 {0{{ajTif)) -. 1.6 {0{{ajT:f)) -^ 1.3 
((9((as/7r)^)). The fact that the g"^ determination of AE falls in increasing orders of 
perturbation theory, whilst the determination from g^ rises, is easily understood from 
the fact that the higher order perturbative QCD corrections to the Bjorken sum rule 
are much larger than those of the singlet combinations of structure functions. The 
main perturbative correction to the non-singlet part comes from the isovector term 
-^Qa which reverses its sign between the neutron and the proton. In the deuteron 
this term is absent and the effect of the remaining og term is small, and therefore 
the Qi determination does not change much in increasing orders of perturbation 
theory. 

We conclude that a consistent overall picture of AE and As is emerging, which 
leads to the values quoted in equation (^ at the beginning of this paper. 

The overall consistency of the different determinations of AE and As that we 
find is not affected by the possible log(l/a;) behaviour |Q of the ^'f'"' that we dis- 



cussed earlier, since this consistency is simply a consequence of the data and the 
assumed correctness of the Bjorken sum rule. However, such behaviour would alter 
the determinations of the individual Ag quoted in equation (|I|). As pointed out 
in [Q , this could shift the experimental values of the integrals Ff '" and hence the 
inferred values of the Ag by about one standard deviation towards a higher value of 
AE and a lower value of As. As also pointed out in [^], these shifts could be much 



larger if the small-x behaviours of g^'^ were even more singular, but this possibility 
does not seem to be very well motivated theoretically. For the reasons discussed 
earlier in connection with the Bjorken sum rule, we have not yet included estimates 
of singlet higher- twist effects 

A (v^^vn^ „ (-0.02 ± 0.01) GeV^ 

in the extraction of the Ag. If we include them in the global fit, we find 

An = 0.85 ±0.03, Arf = -0.41 ± 0.03, As = -0.08 ± 0.03, AE = 0.37 ±0.07 

(12) 
at Q^ = 10 GeV^, and we regard the differences between these and the values in 
equation (p as indicative of the possible theoretical systematic errors. 

We cannot resist commenting that the value of AE that we obtain is much closer 
to the light current-quark model value of zero in the large- iVc and chiral limit [|l^ 
than it is to the most naive quark model of unity, and that finite- A^c and m^ 7^ 



corrections |13],[o0 were estimated to be capable of altering AE by about 0.3 



This picture may be reconciled with otherwise highly-successful models based on 
constituent quarks, if the latter are regarded as effective constructs that contain 
non-trivial internal chiral and gluonic structure [l6li|6^ , |6^ , |6^ , |65| ,|6^ . 



Finally, as an application of these results, we discuss their implications for the 
couplings to nucleons of certain candidates for non-baryonic Dark Matter. We con- 
sider first the lightest supersymmetric particle (LSP), which we assume to be a neu- 
tralino, i.e., some combination of neutral, non-strongly-interacting spin-1/2 partners 
of the SU{2) gauge boson W^, the U{1) gauge boson B, and the neutral Higgs bo- 
son H^2 ill the minimal supersymmetric extension of the Standard Model (MSSM). 
Since the neutrahno has spin 1/2, its couphngs to nucleons include a spin-dependent 
part that is related to the contributions Aq of the different quark species to the nu- 
cleon spin. Experiments at LEP and elsewhere constrain the parameters of the 
MSSM in such a way that the LSP cannot be an approximately pure photino or 
Higgsino. However, it could be an approximately pure U{1) gaugino B, in which 
case its spin-dependent coupling to a proton is proportional to 

ttp = 17/36AM + 5/36(Arf + As) = {6gA + 2a8 + 9AS)/36 (13) 

and the analogous coupling a„ to the neutron is given by a similar formula with 
Au and Ad interchanged, in the limit of a small momentum transfer from the LSP 
to the proton and misp "^ T^q- The couplings ap,„ should be evaluated using the 
value of the AS evolved to a renormalization scale of order uiq, which we take for 
illustration to be 500 GeV. Using the values (and errors) of the Aq given in equation 
(|T]), we find 

ap= 0.32 ±0.017 

(14) 
a„ = -0.10 ±0.017 

We note that the error on Op is only about 5%, whilst the error on a„ is about 17%. 
This implies that one can estimate relatively reliably the spin-dependent coupling of 
the LSP to odd-even nuclei such as ^^K or ^^Nb, whose spins are carried essentially 
by protons, whereas the uncertainty is somewhat larger for even-odd nuclei such as 
^^Ge or ^^Si, whose spins are carried essentially by neutrons. Predictions based on 
the naive quark model [^ for the spin content of the nucleon are very misleading, 



particularly for the neutron. 

Next, we turn to the couplings of the axion to nucleons, which are given by 

Cap = 2 [-2.76 An - 1.13 Ad + 0.89 As - cos 2/3 (Am -Ad- As)], 

(15) 
Can = 2 [-2.76 Ad - 1.13 Am ± 0.89 As - cos 2(3 {Ad - Au - As)] 



for the proton and neutron respectively. In this case, the low-renormahzation scale 
values of the Ag are the relevant ones, and we evaluate them at the scale Q = 1 
GeV, which might be appropriate for the core of a neutron star. We find 

C^p = (-3.9 ± 0.4) - (2.68 ± 0.06) cos 2/5 

(16) 
Can = (0.19 ± 0.4) + (2.35 ± 0.06) cos 2/3 

which can be compared with the naive quark model values given in equation (4) of 



3^ . We see that the axion-proton coupling is relatively well determined, except for 
large values of tan/5, whilst the axion-neutron coupling is more sensitive at interme- 
diate values of tan /?, to the experimental errors in the determination of the Aq. As 
was discussed in |^ and references therein, the dominant axion emission process 



from the core of a neutron star is likely to be axion bremsstrahlung in nucleon- 
nucleon collisions, which is sensitive to a combination of Cap and Can- The relative 
weights of neutron-neutron, neutron-proton and proton-proton bremsstrahlung pro- 
cesses are given roughly by 

CaJ + 0.83{Can + Cap)^ + 0.47^/ (17) 

One tries to bound the axion decay constant /„ by constraining the rate of axion 
emission from the core of the neutron star. Using the numbers given in equation 
(|16D, we find a sensitivity 

Afa/fa = 30% to 50% (18) 

This source of error in the constraint on fa is much less than other sources of error, in 
the nuclear equation of state, for example. We conclude that our present knowledge 
of the Aq is adequate for the purpose of bounding /„. 

We conclude that the available polarized structure function data are in per- 
fect agreement with QCD, once higher-order perturbative corrections are taken into 
account. These are indeed essential: the naive parton model is not adequate to 
describe the nucleon as viewed through polarized lenses. The polarized structure 
function data already yield an interestingly precise value of a^, that future rounds of 
data could render highly competitive with other determinations. The total nucleon 
spin fraction AS carried by quarks is around 30 %, and the strange contribution 
As < 0. The dynamical mechanism that explains these findings remains a puzzle, 
and the challenge to reconcile them with the naive constituent quark model persists. 
More high-precision data would be most welcome, in order to discriminate between 
the chiral soliton and axial U{1) anomaly interpretations of the data, to probe the 
behaviours of (/f " at small x, to explore the Q^- dependence of Ai, and hopefully to 
extract higher-twist effects and AE from data at different values of Q^. Other ex- 
periments, for example on elastic z/iV scattering |Q^,^,0 and using high-energy 
polarized proton beams to measure gluon polarization directly, may also help to 
disentangle the nucleon spin decomposition more completely. This is not only a fas- 
cinating way of probing non-perturbative QCD, but also has implications in other 



areas of physics, for example in searches for Dark Matter particles, as discussed 
above. Floreat nucleon spin physics! 
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Figure Captions: 



Fig. 1. Values of as{Q^ = 2.5 GeV^) extracted from the E142 T^ and E143 T? mea- 
surements using the Bjorken sum rule and increasing orders of QCD perturbation 
theory. The last point includes the estimate of the 4-th order perturbative coefficient 
in 12^ 



Fig. 2. Values of a^iQ^ = 2.5 GeV^) extracted from the E142 F^ and E143 F? 
using the Bjorken sum rule and assuming a given non-zero value of the higher twist 
coefficient Cht (cf. eq. (|^)). Dashed lines denote error bands corresponding to one 
standard deviation with respect to the central value of ag- Vertical dotted lines 
denote the range of Cht given in eq. (0) . 



Fig. 3. The value of asiQ"^ = 2.5 GeV^) from the Bjorken sum rule, using all 
available perturbative QCD corrections but neglecting possible higher-twist effects, 



shown together with a compilation of world data on as{Q) as given in ref. |^ . 



Fig. 4. EMC, SMC and preliminary E143 p data on F?(Q2), together with F?(Q2) 
inferred from the E142 ^He and SMC D data using the Bjorken sum rule. The upper 
continuous curve shows the naively-suggested value holding if As = 0, together 
with an error band plotted with dotted curves. The lower curve shows the prediction 
for F?(g2) that would be obtained if As = -0.10 ± 0.04. 



Fig. 5. The values of AE((5^=10 GeV^) extracted from each experiment plotted as 
functions of the increasing order of QCD perturbation theory used in obtaining AS 
from the data. The last point includes the estimate of the 4-th order perturbative 
coefficient in 
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